Witnessing trustworthy single-photon entanglement 
with local homodyne measurements 
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Single-photon entangled states, i.e. states describing two optical paths sharing a single photon, 
constitute the simplest form of entanglement. Yet they provide a valuable resource in quantum 
information science. Specifically, they lie at the heart of quantum networks, as they can be used 
for quantum teleportation, swapped and purified with linear optics. The main drawback of such 
entanglement is the difficulty in measuring it. Here, we present and experimentally test an entan- 
glement witness allowing one not only to say whether a given state is path-entangled but also that 
entanglement lies in the subspace where the optical paths are each filled with one photon at most, 
i.e. refers to single-photon entanglement. It uses local homodyning only and relies on no assumption 
about the Hilbert space dimension of the measured system. Our work provides a simple and trustful 
method for verifying the proper functioning of future quantum networks. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.50. Dv 



Motivations Quantum networks pQ provide broad ca- 
pabilities, ranging from long distance quantum commu- 
nication at large scales [H [5] , to the simulation of quan- 
tum many-body systems \A- in tabletop implementations. 
Remarkable progresses have been made in practice [1HZ] 
and experimental capabilities are now advancing into a 
domain of rudimentary functionality for quantum nodes 
connected by quantum channels [8-llJ. Surprisingly, the 
task of checking that a newly implemented quantum net- 
work performs well remains non-trivial. 
In the past decade, a great number of architectures based 
on atomic ensembles and linear optics have been pro- 
posed |12) . We now know that quantum networks based 
on single-photon entanglement [13] . i.e. entangled states 
of the form 

-L(|i^|o) B + |o>^|i) B ) (i) 

where A and B are two spatial modes sharing a de- 
localized photon, are very attractive: They require 
significantly fewer resources than the other architec- 
tures and are less sensitive to memory and photon 
detector inefficiencies |12j . Furthermore, they are 
efficient when combined with temporal multiplexing 
[14] . However, such networks have a major drawback: 
The detection of single-photon entangled states is very 
challenging. One cannot resort, for example, to violating 
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a Bell inequality given solely photon counting techniques. 

Hitherto, there are three prescribed methods to detect 
single-photon entanglement. The first one converts two 
copies of a single-photon entangled state into one copy of 
two-particle entanglement. Starting from entanglement 
(|lU|0) fll + |0)^|l) Bl ) ® (\1)a 2 \Q)b 2 + \0)a 2 \1)b 2 ) 
between the modes A\ and B\ and between Ai and B2 , it 
basically consists of a post-selective projection onto the 
subspace with one excitation in each location, yielding 
|1)ai|1)b 2 + |1)a 2 |1)bi EH- The latter is analogous to 
conventional polarization or time-bin entanglement and 
any witness suited for such entanglement can thus be 
used to post-selectively detect single-photon entangle- 
ment. Nevertheless, this approach is not fully satisfying 
conceptually because it relies on post-selection. Further- 
more, for practical implementation, the need to create 
two copies requires twice the number of resources at 
each node. 

The second method is based on partial quantum state 
tomography. Specifically, one reconstructs a reduced 
density matrix that corresponds to a projection of 
the full density matrix into a subspace with at most 
one photon locally. The presence of entanglement is 
then inferred from an entanglement measure computed 
from the reduced density matrix 15J. Specifically, 
this tomographic approach requires the knowledge of 
probabilities p mn of having m photons in mode A and 
n in mode B, where m,n € {0,1}, and the visibility 
V of the single-photon interference pattern obtained 
by combining the modes A & B into a beam-splitter. 
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Although it has triggered highly successful experiments 
[T5HT8] . the approach presented in Ref. [15] cannot be 
directly used in large scale networks when one needs 
to check the entanglement between far away locations, 
since the knowledge of V relies on a joint measurement 
of A & B modes. 

The last method uses local homodyne detections and 
provides, a priori, a full tomography of the state that 
can subsequently be used to measure the entanglement 
[191 120) . However, the tomographic approach requires 
a number of measurements which increases with the 
dimension of the state being measured [2T] . In practice, 
one could be tempted to make an assumption on the 
regularity of the measured Wigner function to reduce 
the number of measurements or, equivalently, on the 
dimension of the system's Hilbert space, especially when 
focusing on single-photon entanglement. But this would 
amount to make an assumption about the system that 
we want to characterize. One can also estimate the 
dimension of the state from measurements but it is not 
clearly established how errors on this estimation can 
affect the conclusion about the presence of entanglement. 
More generally, the exponential increase of required 
measurements with the number of measured subsystems 
makes the tomography not suited to decide on the 
presence of entanglement in quantum networks |22j . 
contrary to entanglement witnesses [25] . 

Principle Here, we propose a simple approach to 
witness single-photon entanglement which relies on lo- 
cal measurements only and needs neither post-selection 
nor assumption on the tested state. The basic prin- 
ciple is drawn in Fig. [T] Two distant observers, Al- 
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FIG. 1: Principle of the proposed entanglement witness. 

ice & Bob, share a quantum state. To check whether 
it is entangled, each of them randomly chooses a mea- 
surement among two quadratures, {X, P} for Alice and 
{X + P. X — P} for Bob. At each run, they obtain a 
real number. They then process the results to get binary 
outcomes using a sign binning, i.e. they attribute the 
result —1 if the result is negative and +1 otherwise. By 
repeating the experiment several times, Alice & Bob can 
compute the conditional probabilities p(a,b\x,y) where 
a, b G {-1, +1}, x G {X, P} and y G {X + P, X - P}. 
Substituting these probabilities by their values into the 
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P,X-P, 



(2) 



where E XtV = Y, a ,b={-i,+i} P( a = b \ x >v) - P( a ^ b\x,y), 
they obtain a real number S. 

The value of S can easily be obtained under the as- 
sumption that Alice & Bob have each a qubit. Indeed, 
in the Fock basis {|0),|1)}, the measurement of the 
X quadrature with sign binning is equivalent to a 
noisy a x measurement [25l[26], j'^di \x)(x\ — / °° dx- 



\x)(x\ = y and similarly, the P quadrature 
corresponds to o~ y with the same noise. For the setting 
choice {X, P} and {X + P,X - P}, the state Q thus 
yields to S = 2y2x - ps 1.8. Furthermore, the maximum 
value that can be obtained with a separable state belong- 
ing to the subspace {|0), \l)}® 2 is S scp = V2 x f « 0.9 
[27]. Since S is smaller than 2, the proposed CHSH-likc 
test does not highlight the non-local characteristic of 
a single photon delocalized among two modes, but it 
does provide an attractive entanglement witness: If the 
measured CHSH value is larger than S sep , Alice & Bob 
can conclude that they share an entangled state. 
This holds for qubits only. In practice, however, the 
state describing the modes A & B includes multi-photon 
components, and does not reduce to a two-qubit state. 
We show below how the entanglement witness can be 
extended to the case of arbitrary dimensional bipartite 
states. First, we show how Alice & Bob can accurately 
estimate the probability that their state lies out of a 
two-qubit space {|0),|1)}® 2 . We then demonstrate that 
this probability can be used to upper bound the maximal 
CHSH value that can be obtained with separable states. 

Bounding the Hilbert space dimension Let us con- 
sider the case where Alice & Bob do not have qubits, 
but quantum states of arbitrary dimension. First, they 
need to bound the probability that at least one of their 
modes is populated with more than one photon p(riA > 
2 U ng > 2). This can be realized without assumption 
on the Hilbert space dimension, by first determining the 
probabilities p(riA = j) {p{nB = j)) of having j pho- 
tons in Alice's (Bob's) mode using local homodyning with 
phase averaged local oscillators through a direct integra- 
tion of the obtained data with a pattern function |28) . 
The joint probability p{nj± > 2Un# > 2) = p(ua > 
2) + p(ns > 2) — p{ua > 2 n ub > 2) can then be 
bounded by the parameter p* defined as follows 



,p(nA = j)+p{n B = j) 



(3) 



We now show how the knowledge of p* can be used 
to construct an operational witness for single-photon 
entanglement. 
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FIG. 2: Separable bound S^p* as a function of the probability 
that at least one of the two protagonists gets more than one 
photon pfjlA > 2 U hb > 2). The blue curve allows one to 
know the maximum value of the CHSH polynomial S™p* that 
a state separable in the {|0},|1}}® 2 subspace can reach. If 



Sobs ^ S*. 



one can conclude that the projection of the state 



in the subspace with zero and one photon locally is entangled. 



Evaluating the maximal CHSH value with separable 
states Consider the general case where p* ^ 0, i.e. 
p(n-A > 2 U us > 2) ^ a priori. The state of Alice 
& Bob can be described by the density matrix 



Pn, A <\C\n B <\ Pcoh 

^coh Pn A >2Un B >2 



(4) 



where p nA <mn B <i denotes the 4x4 block with at most 
one photon per mode, p nA >2un B >2 refers to the block 
where at least one of the two modes contains at least two 
photons and p co h is associated to the coherence between 
these two blocks. Since p nA >2un B >2 possibly spans a 
Hilbert space of infinite dimension, there could be an in- 
finite number of coherence terms. However, a few of them 
give a non-zero contribution to the CHSH polynomial if 
a phase averaged homodyne detection is used at each lo- 
cation. Specifically, consider the case where Alice & Bob 
perform the measurements X VA = cos (paX + sin^^P 
and X, 



<PB 



cos tpsX + sin (psP respectively, where ipA 
and (fiB are random variables such that (e zklpA - B ) = 0, 
k £ N* but where the phase difference <pa ~ <Pb = A(p is 
fixed. This only requires classical but not quantum com- 
munication, and hence can only decrease the entangle- 
ment that Alice & Bob potentially share. In particular, 
if Alice can choose a measurement among the two quada- 

tures \X, n i , X,„2 ) and if Bob's choice reduces to one of 

r A v A J 

the quadratures {X^, X^} such that ip A — ip B = — j , 

"Pa ~ <P% = f ><Pa ~ <Pb = f and "Pa ~ 9% = x> we show 
in the appendix that the CHSH polynomial correspond- 



ing to the state Q is bounded by 
16 
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(5) 



- (^[(20|p coh |ll) 
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(02|p coh |ll) 



-2V2p(n A >2Un B > 2) 



where 5ft denotes the real part. For a given 
value of p{n A > 2 U n B > 2), S m& * can be di- 
rectly maximized over the set of physical states 
(p e {|0),|l),|2)}® 2 ,tr(p) < l,p > 0) that sat- 
isfy the observed photon number distributions, i.e. 
Poo + Pio + Poi + Pn = 1 - p(n>A > 2 U n B > 2) and that 
are separable in the {|0), |1)}® 2 subspace, i.e. for which 
the projection into this subspace remains positive under 
partial transposition (PPT) [HI [3D]. Figure [2] shows the 
result of this optimization 5™p X . 

Witnessing single-photon entanglement This pro- 
vides a truly state-independent witness [3T] of en- 
tanglement: First, the protagonists determine the 
local photon-number distributions, from which they 
deduce an upper bound p* on the joint probability 
p(n A > 2Un B > 2). Secondly, they deduce S™ x {p*) 
form Fig. [2] (See the appendix). Thirdly, they measure 
the CHSH value S b s by randomly choosing measure- 
ments among {X Va , X^ } and {X^^X^} respectively 
and by subsequently computing the CHSH polynomial 
through Eq. If S ohs > S™ x (p*), Alice & Bob 

know that the projection of their state into the subspace 
{JO),!!)}® 2 has a negative partial transpose, i.e. they 
can safely conclude that the state is entangled and that 
the entanglement resides in the subspace with at most 
one photon locally. 

Importantly, a tighter bound can be obtained if 
gmax j g maximized over the set of states with 
a positive partial transpose not only satisfying 
Poo + Pio + Poi + Pn = 1 - p(riA > 2 U n B > 2) 
but also reproducing the locally measured probabilities 
p(nA = j) (p(riB = j)) for having j = 0, 1 photon in 
Alice's (Bob's) mode. These additional constraints have 
been taken into account for the computation of the 
separable bounds related to the experiment presented 
below (see appendix). 

Proof- of-principle experiment We start off with a her- 
alded single photon generated by a conditional prepa- 
ration technique operated on a two-mode squeezed vac- 
uum emitted by a type-II optical parametric oscillator 
|32) . Without correction for detection loss, the overall 
fidelity reaches 70%. Single-photon entanglement is ob- 
tained by sending the created photon into a beam-splitter 
[33] . Specifically, by controlling the angle 9 of a half- wave 
plate relative to the axis of a polarizing beam-splitter 
(PBS), we create a tunable single-photon entangled state 
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cos(20)|O)a|1).b + sin(26»)|l) A |0) B between the two out- 
put modes of the PBS, as sketched on Fig. [3] 
The local oscillators that Alice & Bob need to reveal 





FIG. 3: Experimental setup. A tunable single-photon entan- 
gled state is created by sending a heralded single photon on 
a tunable beam-splitter based on a polarizing beam-splitter 
(PBS) and a half-wave plate (A/2). The proposed witness is 
then tested with two independent homodyne detections (Alice 
& Bob). The local oscillator is superposed to each modes via 
the first PBS. Its global phase is swept with a piezoelectric 
actuator. The relative phase Aip is set with a combination of 
birefringent plates. 



entanglement are obtained by impinging a bright beam 
on the second input of the PBS: one polarization mode 
is used as the quantum channel, the orthogonal one con- 
veys the local oscillator, ft is worth mentioning that this 
technique can be implemented over long distances as real- 
ized in field implementation of quantum key distribution 
[3~4"j . The relative phase between Alice's and Bob's local 
oscillators Atp is fixed by choosing an appropriate ellip- 
tical polarization of the bright beam just before the PBS 
[35] . In practice, the setting difference is calibrated by 
observing the dephasing of interference fringes (the quan- 
tum state is replaced here by a coherent state) . A global 
phase averaging is obtained by sweeping a piezoelectric 
transducer located on the path of the bright beam before 
the PBS. 

For each heralding event, Alice & Bob each obtain a real 
valued outcome which is extracted from homodyne pho- 
tocurrents. Accumulating 200000 events for each quadra- 
ture relative measurements, they deduce the value of the 
CHSH polynomial Sobs- The same homodyne measure- 
ments also provide the local photon number distributions 
which are used to compute the separable bound S™*. 
We emphasize that the separable bound is here obtained 
by maximizing the CHSH value over the set of separable 
states that fulfill the locally measured photon-number 
occupation probabilities p(ua — j) and p{ng — j) for 
having j = 0, 1 photon in Alice & Bob's mode respec- 
tively. Furthermore, it takes several errors into account, 
for example, errors related to the quadrature measure- 



ment imperfections were considered (see appendix). The 
procedure is repeated for various angles 9 ranging from 
to 45°. Fig. [4] shows the main result, i.e. the ob- 
served CHSH values and the separable bounds as a func- 
tion of 9. One sees that they both reach maximal values 
around 9 = 22.5° where Alice & Bob ideally share a max- 
imally entangled state. The small deviation between the 
observed value S ohs (9 = 22.5°) « 1.33 and the CHSH 
value that would be obtained with a maximally entan- 
gled state (1.8) demonstrates that the overall source and 
detection efficiencies are very high. Furthermore, the ob- 
served CHSH values are almost all larger than the separa- 
ble bounds when dealing with entangled states (9 ^ 0°,- 
45°). This shows the great robustness of the proposed 
witness. 
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FIG. 4: Observed CHSH values Sobs (the size of points ac- 



scp 

> S™ x , one 



counts for statistical errors) and separable bound S, 
function of beam-splitter angle 6. When S b s > S L 
can conclude that the measured state is entangled. 

Conclusion We have presented and experimentally 
tested a witness for single-photon entanglement which 
does not need postselection, uses local measurements 
only and does not rely on assumptions about the dimen- 
sion of the measured state. Note that our witness can 
be easily adapted to detect few-photon entanglement 
without additional complications. We believe that it 
will naturally find applications in long distance quantum 
communication, allowing users to check whether two 
remote nodes of a given quantum network are entangled. 
One important challenge in this context is to reveal 
the entanglement shared by a large number of parties. 
Finding Bell inequalities that could be used as witnesses 
for multi-partite single-photon entanglement is work for 
future. 
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Apendix A: Witnessing single-photon entanglement 
in qudit spaces: Theory 



As mentioned in the main text, if Alice & Bob can 
guarantee that they perform measurements on qubits, 
they can demonstrate entanglement by simply perform- 
ing a Bell-like test by randomly choosing a setting among 
two quadratures, {X, P} for Alice and {X+P, X — P} for 
Bob. If the resulting CHSH value is larger than 0.9, they 
can safely conclude that their state is entangled. In prac- 
tice, however, it is challenging to show that the systems 
one is measuring are well described by qubits. Below, we 
detail the procedure to follow in the general case where 
Alice & Bob have qudits. 



CHSH polynomial for phase averaged homodyning 



Let 



^2 c ijki\ij)(kl\ 

ijkl 



be the state that Alice & Bob share, \ij) describing the 
state with i photons in Alice's mode and j photons in 
Bob's one. Consider that the phase of local oscillators 
that are required for homodyning, is averaged, i.e. Alice 
& Bob perform the measurements 

x v t A = cos Va x + sin Va p , 
X^i — cos (figX + sin ip e B P 

respectively, with ip A and <pg random variables satisfy- 

■ i ■ i 

ing (e lVA ) — (e %VB ) — 0. Further consider that the phase 
difference Aip u = ip l A — ip l B to be tunable such that Al- 
ice can choose a measurement among the two quadatures 
{X^, X^} relative to Bob's choice {X v i g , }. (This 
can only underestimate the entanglement in average be- 
cause this can be realized by local operations and clas- 
sical communications. Namely, they could each apply 
a random phase shift to a shared local oscillator such 
that the phase difference is fixed, subsequently choose 
a quadrature measurement locally and ignore the infor- 
mation about the individual phase shifts). At each run, 
Alice & Bob obtain each a real number. They then pro- 
cess the results to get binary outcomes using a sign bin- 
ning, i.e. they assign the result —1 if the real number is 
negative and +1 otherwise. By repeating this procedure 
many times, they can access the probability that both 
get +1 for instance, knowing that they chose X^i and 



p(+l,+l|A> A ,A> s ) =£( e V B (^-(fc+*)) )c 



ijkl 



ijkl 



xe 



iA v il (i-k) 



/•oo 

/ dx4>i(x)4>k(x) 
Jo 



x / dy<i>j{y)(t>i{y) 

Jo 

where (f>i(x) = (x\i). One sees that the off-diagonal el- 
ements \ij)(kl\ with different numbers of photons (i + 
j 7^ k + I) do not contribute to the probabilities 
p(a,b\X ip i ,X i), with a, b = {— 1,+1}. Furthermore, 
since for all n and m having the same parity, 

/•oo A) 

/ dy<j) n {y)4> m {y) = j dy4> n {y)<j> m {y) = -<5„, m 

JO J-oo 1 

the terms Cijki for which either i-k or j-l is an odd num- 
ber, have a zero contribution to the correlator 

E x^x^ =E., H -i,+i} P( a = b \ X ^ X v %) 

- V {a + b\X^X^ B ). 

Finally, one easily checks that the remaining c^ki terms 
satisfying 



i + j = k + I, 

\i — k\ = 1 mod 2, 

\j -l\ = l mod 2, 



(6) 
(7) 
(8) 



yield 



E 



C-ijkl 



ijkl 



/•oo 

/ dxct)i{x)4) k (x)e 
Jo 

/•OO 

/ dy<f>j(y)<f>i(y). 
Jo 



Specifically, if Atp 11 = -f , Aip 12 = \ , Atp 21 = f and 
Aip 22 = ^f, the value of the CHSH polynomial 

S = Ex j ,x i + Ex ! .x 2 + Ex 2 ,x i — Ex 2 ,x 2 

fA <?B f>A *B fA *B fA *>B 

is given by 

S= 16 Cjjki cos ((i - k) 

ijkl 

/•oo />oo 

x / dx<j>i{x)<j> k {x) . I dy<f>j(y)(j)i(y). (9) 
Jo Jo 



Optimizing the CHSH value over the set of 
separable states 

Without loss of generality, Alice & Bob's state can be 
written as 



. Pn A <lnn B <l Pcoh 

P = I "t 

Pcoh Pn A >2Un B >2 



(10) 
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where p nA <mn B <i denotes the block with at most one 
photon per mode, p nA >2un B >2 refers to the block where 
at least one of the two modes is populated with more 
than one photon and p co h is associated to the coherence 
between these two blocks. Taking the constraints (|6]),([7]) 
and Q into account and using the formula the cor- 
responding CHSH polynomial reduces to 



16 
" V2tt 

7T V 



(01 1 A 



n A<lC\n & < 1 



(20| A 



'coh 



11) 



|10) 
"(02|p, 



(11) 



coli 



11) 



where 3? denotes the real part and S Pn >2UnB>2 is associ- 
ated to the CHSH value obtained from p nA >2Un B >2- The 
goal is now to optimize S over the set of separable states. 
Let p(n A > 2 U n B > 2) be the probability that at least 
one of the two protogonists has more than one photon, 
i.e. 



tr Pn A >2un B >2 = p(n A > 2 U n B > 2). 



(12) 



Since p nA >2un B >2 cannot be obtained by local pro- 
jections, it may maximally contribute to CHSH, i.e 
S. 



,,„ , _, 2 < 2^2 x p(n A > 2 U n B > 2). S is thus 

"ed by 

(13) 



upper bounded by 
S < 5 max 



(p(n A > 2Un B > 2)) 

16 r 

(01|/?n A <lnnB<l| 10 ) 



-5ft 



V2ir 
+ - (jft[(20|p coh 

7T \ L 

+2y/2 x p(n A > 2 U n B > 2). 



11) 



ift 



(02|p coh |ll) 



Needless to say, separable states are physical states. 
They are thus represented by positive matrices with 
a trace (tr) equal to one. Furthermore, the Peres- 
Horodecki criterion [TJ |5] states that for any separable 
state /9 sop , its partial transpose p^? p has non- negative 
eigenvalues. The optimization of <> max over the set of 
separable states in the subspace with at most one pho- 
ton locally is thus a problem that can be summarized as 
follows 



max 
P e{|o>,|i>,|2>}®2 

s. t. 



S max (p{n A > 2Un B > 2)) 

P > 
tr(p) < 1 



(14) 



IMI 

v.0/1 0/1, 
1 

i,j=0 



> 



1 -p(n A >2Un B >2). 



The constraint tr(p) < 1 comes from the fact that the 
optimization is performed over finite dimension (9 x 9) 



matrices that can either represent a physical state or 
that can be obtained by local projections of states 
spanning Hilbert spaces with a larger dimension. Jlo/i 
stands for the projection into the subspace with at most 
one photon locally. Eq. ( |14[ ) is a linear optimization with 
scmidcfinite positive constraints which can be efficiently 
solved numerically [3]- The result of the optimization 
*Ssep X i s shown in Fig. 2 (main text) as a function of 
p( n A > 2 U n B > 2). It is also given in Fig. [5] (dashed 
line). If a physical state that satisfies the last condition 
leads to a S value larger than S™ ax , one concludes that 
the partial transpose has at least one negative eigenvalue 
and hence, the state is entangled and entanglement lies 
in the subspace with zero and one photon locally, i.e. 
refers to single-photon entanglement. 




0.02 0.04 0.06 o.c 

p(n A >2Un B >2) 



FIG. 5: Separable bound S^p* as a function of the probability 
that at least one of the two protagonists gets more than one 
photon p(n A > 2Unj > 2). The dashed curve allows one 
to know the maximum value of the CHSH polynomial S^p"* 
that a state that is separable in the {|0), |1}}® 2 subspace can 
reach. If the observed CHSH value Sobs > Sscp*, one can 
conclude that the projection of the state in the subspace with 
zero and one photon locally is entangled. The full curve give 
the maximum value that a separable state can reach. If the 
the oberved CHSH value is higher than the latter, one can 
conclude about entanglement but without saying where the 
entanglement lies. 



It is worth mentioning that if the observed CHSH 
value is not higher than the separable bound obtained 
through the optimization presented in Eqs. (14 1, one 



can still conclude about the presence of entanglement. 
Indeed, the optimization can be performed over the 
set of separable states, but not only the ones that are 
separable in the subspace with at most one photon 

locally. The constrain (ilo/i ^Ilo/i) — simply 
needs to be replaced by p Tb > 0. This leads to the second 
separable bound (full curve) presented in Fig. [5] If the 
observed CHSH value is larger than the bound obtained 
in this manner, one can conclude about the presence 
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of entanglement but we cannot say in which subspace 
entanglement lies. 



Note also that the optimization ( 14 ) does not take the 
individually measured probabilities pin a — j) (p(«s = 
j)) for having j photons in Alice's (Bob's) mode into 
account. These probabilities add constraints reducing 
the size of the set of separable states over which g max 
is optimized and thus, provide a tighter bound S™*. 
We use them for the calculation of the separable bounds 



related to the experiment (see the optimization (19)) 



General procedure to follow for witnessing 
single-photon entanglement 

Below, we present the procedure to follow in order to 
conclude about the presence of entanglement with the 
proposed witness. 

• Firstly, the local photon number distributions 
p{n a = 0),p(n b = 0),p(n a = l),p(n b = 1) are 
accessed using local phase-averaged quantum state 
tomography. 

• The joint probability p(n a > 2(Jn b > 2) defined by 
p(n a > 2) + p(n b > 2) - p{n a > 2 n n b > 2) 

is then upper bounded by 

P* = 2 - p(n a = 0) - p(n a = 1) 
-p(n b = 0) -p(n b = 1). 

• Thirdly, p* is used to obtain the separable bound 
5™p X (p*) deduced from the following optimization 

max : S roajc (p*) 
P e{|0>,|i>,|2>}®2 

s. t. : p > 

tr(p) < 1 

IMI >o 

.0/1 0/1/ 



1 

E 



Pij >l-p*. (15) 



Alternatively, S™p X (p*) can simply be obtained 
from Fig. 2 (main text) because S™p X is a mono- 
tonically increasing function of p(n a > 2 U n b > 
2). Both methods provide an upper bound on 
S™ x (p(n a > 2 U n b > 2)) < S™*(p*). 

• Fourthly, the CHSH value is measured. In princi- 
ple, this is done by proposing Alice and Bob to ran- 
domly choose measurements among { X^i , X^i } 



and {X^i , X^ } respectively and to subsequently 
compute the CHSH polynomial. However, since 
Ex i ,x i = —Ex 2 ,x 2 an d Ex t .x 2 = 
Ex 2 ,x i j &abs can be obtained in practice from 
the measurements of two correlators only 



"Sob 



• Finally, 
If S obs 



2E 



x 



,x 



2E 



x 



(16) 



Sobs 
> 



and 

crnax 
^sep 



Cm ax 
^sep 



(p*) 



(p*) are compared, 
(which implies that 



S obs > S™*(p(n a > 2Un b > 2))), one con- 
cludes that the measured state is entangled at the 
single-photon level. Otherwise, we cannot form a 
conclusion, the state can either be separable or 
entangled. 



Appendix B: Experimental details 
Heralded creation of single-photon entanglement 

A continuous-wave frequency-doubled Nd:YAG laser 
(Diabolo, Innolight) pumps a triply-resonant type-II 
phase-matched optical parametric amplifier based on a 
KTP crystal to generate below threshold, a two-mode 
squeezed state [B]. The created modes are orthogonally 
polarized and are deterministically separated at the out- 
put of the amplifier. One of the two modes is send to 
a single-photon detector (superconducting single-photon 
detector, with a quantum efficiency of 7% at 1064nm) 
after filtering of the non-degenerate modes. A detection 
event ideally heralds the generation of a single-photon 
state on the twin mode. The heralding rate is around 
30 kHz. (See [5] for more details about the source). By 
controlling the polarization of the heralded photon via a 
half- wave plate and by subsequently sending it into a po- 
larizing beam splitter (PBS), one gets a versatile source 
producing states of the form 

\iI>(0))ab - cos(20)|O)a|1)b + sin(20)|l) A |O) B (17) 

where 9 is the angle of the half-wave plate relative to 
the axis of the PBS. For 9 = 0°, the created state is 
separable, whereas for 9 = 22.5° it becomes maximally 
entangled. This source is thus particularly well suited to 
test the proposed entanglement witness. 



Homodyning 

Each spatial mode is then detected using an indepen- 
dent homodyne detection. The required local oscillators 
are obtained by impinging a bright beam into the second 
input of the PBS. The phase of Alice's local oscillator 
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is controlled relative to Bob's one by choosing appro- 
priately the polarization of the bright beam just before 
the PBS [6]. To meet the witness requirements, a phase 
averaging is realized by sweeping a piezoelectric trans- 
ducer located on the path of the bright beam before the 
PBS. The overall efficiency of each homodyne detection 
is 85%, including the quantum efficiency of the photodi- 
odes (Fermionics 500), mode overlap and electronic noise 



Data acquisition 

For each heralding event, Alice and Bob perform a 
quadrature measurement and the corresponding result is 
extracted from homodyne photocurrents. 200000 events 
are accumulated for each relative phase Aiyj^g = ~ and 
A^g = — f • S bs is then deduced from Eq. (16). The 
same results are also used to compute the local photon 
number distributions by phase-averaged quantum state 
tomography using the method given in [8]. We remind 
that this method does not require assumptions on the 
the dimension of the measured state. Table [6] shows the 
results for various angles 0. Alternatively, a MaxLike al- 
gorithm [5] can be used to access the local photon number 
distribution, albeit with a truncation of the Fock space. 



random set of 200 data points for each local probability 
p(n A = 0),p(n B = 0),p(n A = l),p{n B = 1) . . . and the 
corresponding standard deviation provides the desired 
error Sp(n A = 0), Sp(n B — 0), Sp(n A = 1), Sp(n B = 1) . . . 



We now show how the errors on the measurement an- 
gles e and on the local photon-number probability dis- 
tributions 6p(riA = i),Sp(n B — j) have been taken into 
account in the calculation of the separable bound 5™p X . 
Since, in practice, it is necessary to measure two corre- 
lators only, the separable bound can be calculated under 
the assumption that the measurements that have been 
performed are such that 



,11 



4 



3tt 

T 



21 

AB 



7T 

4 



,12 



i.e. that the errors are the same for A^^g and A<p A 2 B 
and for A(p\ 2 B and /S.Lp 2 AB . These erroneous measurement 
angles yield to an upper bound <> max (which replaces the 



one in Eq. (13)) given by 



Error estimations 



Several kind of errors should be taken into account. 

Firstly, statistical errors affect the measured value of 
the CHSH polynomial S &s- These errors are estimated 
in a standard way by using the central limit theorem. 
They are basically very small (see Table |6| because they 
were deduced from 200000 results. 

Secondly, the accuracy with which the relative phase 
Aip AB between Alice & Bob's measurements is estimated 
to be 2°. This means that when Alice and Bob choose 
measurement settings, e.g. corresponding ideally to the 
quadratures X^-X^i , the relative phase A(p AB is not 
exactly equal to \ as it should be, but Aip AB = | + e 11 
where —1° < e 11 < 1°. This error is taken into account 
into the separability bound. 

Thirdly, errors also affect the local photon-number 
probability distributions, which are estimated using the 
phase-averaged homodyne measurements [8] . Evaluating 
these errors is a not trivial task. We use the following 
method for estimating the overall error. The tomog- 
raphy that we use to access the local photon number 
distributions leaves us with the diagonal elements of an 
estimated density matrix /9 os timate- Using this matrix, 
we simulate the quadrature data, then reconstruct 
the diagonal elements of a simulated density matrix 
Psimui- The approach is repeated 200 times, always 
from the same initial state p e stimate- This generates a 



S < S niax (p(n A >2Un B > 2), e 11 , e 12 ) 

= - fsR[(10|p|01)l x (7 - 3? [<10|p|01> 
7r V L J L 

4 ~[(2G|p|ll) + (ll|p|Q2)l x C 



V2tt 

-2V2 x p(n A > 2 U n B > 2) 



(20|H11) + <11H02) 



x D 



where 9 denotes the imaginary part and 



C = 2(cos(e n --)+cos(e 12 
D = 2 (sin(e n - |) + sin(e 12 - 



SO 



x D 



(18) 



gmax can now jj e pti m j zec i over the set of separable 
states, as before. If one takes into account the local 
photon number probability distribution and the errors as 
well, one ends up with the following optimization prob- 
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FIG. 6: Results of local measurements versus the angle of the half-wave plate relative to the axis of the PBS. The table on the 
left gives the photon number distribution in Alice's location. (p_0, p_l, p > 1 denote the probability that zero p(riA = 0), one 
P(tia = 1) and more than one photon p(rhA > 2) occupy Alice's mode. The second table gives Bob's results. The two last tables 
give an upper bound on the probability that at least one of the two protagonists gets more than one photon p(tia > 2Un_g > 2) 
and the observed CHSH value Sobs-) 



lcm 

max : S max (p\ e 11 , e 12 ) (19) 
s.t. : p > 

tr{p) < 1 

imi) 

0/1 0/1/ 

-1° < e n ,e 12 < +1° 
Poo + Pio + Pai +Pn > 1 ~P* - Sp* 
Pqo + Poi + P02 < p(n A = 0) + Sp(n A = 0) 
Pio + Pn + P12 < p(n A = 1) + Sp(n A = 1) 
P20 + P21 + P22 < p(n A > 1) + 5p(n A > 1) 
Poo + Pio + P20 < p{n B = 0) + Sp(n B = 0) 
Poi + Pn + P21 < p{n B = 1) + Sp(n B = 1) 
P02 + P12 + P22 < p{n B > 1) + Sp(n B > 1) 

In addition to errors, we emphasize that the previous 
optimization uses the knowledge of each local photon 
number probability p(n A = 0), p(n A = 1)... This 
provides a tighter separable bound as compared to the 
optimization summarized in Eq. ( |15[ ). 
Note that the optimization ( |l9[ ) is non-linear in e 11 and 
ei2- Hence, the result cannot be obtained with standard 
scmidcfinitc solvers. However, one can convince oneself 



that the maximum value of 5™p X (p*, e 11 , e 12 ) is obtained 
for the extremal choice e 11 = 1°, e 12 = —1°. This 
can be understood intuitively since this choice brings 
the settings closer to each other and thus, helps in 
increasing the CHSH value. Note that the linearity of 
the optimization is recovered once the values of {e , e 12 } 
are fixed. The results of the optimization for e 11 = 1° 
and e 12 = — 1° are shown in Fig. 4 (main text). 
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